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The theoretical understanding of the spin-orbit coupling (SOC) effects at LaAlOs /SrTiOs inter- 
faces and SrTiOs surfaces is still in its infancy. We perform first-principles density-functional-theory 
calculations and derive from these a simple tight-binding Hamiltonian, through a Wannier function 
projection and group theoretical analysis. We find striking differences to the standard Rashba the- 
ory for spin-orbit coupling in semiconductor heterostructures due to multi-orbital effects: by far the 
biggest SOC effect is at the crossing point of the xy and yz (or zx) orbitals; and around the F point 
a Rashba spin splitting with a cubic dependence on the wave vector k is possible. 

PACS numbers: 73.20.-r, 73.21.-b, 79.60.Jv 



Introduction: Very recently, Caviglia et aL[l] and Ben 
Shalom et aL[2] studied the magnetotransport proper- 
ties of the high mobility two-dimensional electron gas 
(2DEG) at the interface between two insulating per- 
ovskite oxides LaAlOg (LAO) and SrTiOs (ST0)[3]. 
They found strong spin-orbit coupling (SOC) effects, 
whose magnitudes can even be tuned by gate voltages. 
This strong SOC has been the basis for many subsequent 
theoretical and experimental studies [4-9]. Despite its in- 
creasing importance, a clear physical picture of the SOC 
effects in LAO/STO interfaces is still missing. 

SOC generally originates from the relativistic correc- 
tion {h/2mlc^){\/V X p) ■ sto the Schrodinger equation, 
with rrie being the free electron mass and V the poten- 
tial in which the electrons move with momentum p and 
spin s. If V{r) has spherical symmetry like in atoms 
or approximately in solids, the form can be reduced to 
^(r)/ • s, where ^(r) denotes the strength of the atomic 
SOC, and / and s are the orbital and spin angular mo- 
menta of the electron. Thus SOC lifts orbital and spin 
degeneracies. For a cubic perovskite such as STO, the 
six initially degenerate t2g orbitals at the T point are in- 
deed split by SOC [10-13], but Kramers degeneracy (time 
reversal symmetry) is preserved. Because bulk STO has 
both crystal inversion symmetry £(/c,t) = £(— /c,t) ^ind 
time reversal symmetry £:(^, t) = the energy e 

at wave vector k is still spin-degenerate. 

This changes for a LAO/STO interface since the inver- 
sion symmetry is broken so that the SOC lifts the spin 
degeneracy of the 2DEG. This two dimensional SOC ef- 
fect is known as Rashba spin splitting, which has been 
widely studied in semiconductor heterostructures [14, 15] 
and metal surfaces [16]. Assuming a nearly free 2DEG 
with effective mass m, its two spin components will 
be split by Ar = 2a nk^ where the Rashba coefficient 
c^R = {h / Am? c^)dV {z) / dz depends on the potential gra- 
dient in the z direction (perpendicular to the interface). 
Naively, since an electric field gives rise to an electrostatic 
potential gradient, A^^ seems to be simply proportional 
to the gate voltage. However, a typical electric field in 



experiments is ~100 V/mm yielding according to the for- 
mula above A^^ 10~^meV[14], which is much smaller 
than the measured values of the order of meV[l, 2]. This 
is because the assumption of nearly free 2DEG, which ig- 
nores the region of ion cores and the asymmetric feature 
of the interface wavefunction, is too simple. In reality, the 
expression for A^^ is much more complicated [17-21] and 
usually treated as a fitting parameter in semiconductor 
heterostructures and metal surfaces. 

The aim of this Letter is the theoretical description 
of SOC in oxide heterostructures and surfaces. From 
density-functional-theory (DFT) calculations we derive 
a tight-binding (TB) Hamiltonian for the low energy t2g- 
orbitals and their spin-splitting. In particular, we show 
that besides the standard /c-linear Rashba spin splitting, 
there can be a /c-cubic spin splitting around the T point 
due to multi-orbital effects. A much larger spin splitting 
occurs at the crossing point of the xy and yz (or zx) 
orbitals. 

Method: Besides bulk STO, we calculate (i) LAO/STO 
(1.5/6.5 layers) [22, 23] which is symmetric with two n- 
type interfaces, (ii) LAO/STO (4/4) and (1/1) which 
is asymmetric with n- and p-type interfaces [24], (iii) 
vacuum/LAO/STO (3/1/4) which has a single n-type 
interface[25, 26], (iv) vacuum/STO (3/7.5) with a SrO 
terminated surface. We fix the in-plane lattice constant 
of the supercells to the calculated equilibrium value of 
STO, and optimize the internal coordinates. The DFT 
calculations have been done using the WIEN2k code 
with the generalized gradient approximation (GGA)[27, 
28]. The SOC is included as a perturbation using the 
scalar-relativistic eigenfunctions of the valence states. 
Through a projection onto maximally localized Wannier 
orbitals [29-31], we construct a realistic TB model, and 
in particular, we develop a way to describe the interface 
asymmetry. 

DFT results for STO: Bulk SrTiOs is a band insula- 
tor with an energy gap between occupied 02p bands and 
unoccupied Tis^^ t2g {yz, zx, xy) bands. The t2g band 
structure calculated by DFT is shown in Fig. 1(a). In 
the absence of SOC, the three t2g orbitals are degenerate 
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FIG. 1: Band structure of i2g orbitals in bulk SrTiOs cal- 
culated by (a) DFT and by (b) a TB model derived in this 
paper. In the absence of spin-orbit coupling, yz^ zx and xy are 
degenerate at the F point. SOC splits the six- fold degenerate 
orbitals into F^ and F^ states separated by Ao =29 meV. 



at F(0, 0, 0) due to an octahedral Oh crystal field around 
the Ti atoms in a perfect perovskite structure. Around F 
the bands can be fit by a parabolic function of the form 
fi^/c^/2m* where the effective mass m* depends on the 
orbital and the direction. 

Along the F— X(7r,0,0) direction (here in units of 1/a 
with a=3.92A being the calculated lattice constant of 
STO), the yz band has a small energy dispersion corre- 
sponding to a heavy mass of 6.8me. In contrast, the zx 
and xy orbitals have the same, large energy dispersion 
and a light effective mass of 0.41me. Including the SOC, 
the six-fold degenerate orbitals are split into a doubly 
and a four-fold degenerate level with an orbital splitting 
of Ao = 29 meV [13], see Fig.l. The energy dispersion of 
the resulting orbitals is now considerably different from 
the initial orbitals. Consequently, the corresponding ef- 
fective masses around F are changed to 1.39, 0.41 and 
0.53me, respectively. 

TB Hamiltonian for STO: To understand the DFT re- 
sults, we use a Wannier projection to obtain the local 
energy and hopping terms of the t2g orbitals. With- 
out SOC, the constructed three-band TB Hamiltonian 
of bulk STO i^o can be expressed in the t^g basis in a 
matrix form: one of the diagonal terms is {xy\HQ\xy) = 
— 2ticoskx — 2ticosky — 2t2Coskz — 4:tscoskxCosky] the 
other two follow by exchanging the x^y^z indices. The 
local energy term is Sq = 3.31eV for all three orbitals. 
The large hopping ti = 0.277eV stems from the large xy 
intra-orbital hopping integral along the x and y direction; 
it is due to two lobes of xy orbitals at the nearest neighbor 
sites pointing to each other along the two directions. In 
contrast, ^2 = 0.031eV and = 0.076eV indicate a much 
smaller hopping integral along the z and (1,1,0) direc- 
tion respectively. We find all orbital-off-diagonal (inter- 
orbital hopping) terms to be negligible. The TB energy 
dispersion for the three orbitals is plotted in Fig. 1(b), 
showing a good agreement with the DFT results. 

We include SOC at the atomic level as = 
^/ • s, where ^ is the atomic SOC strength and de- 
pends on atomic numbers. Under the Oh crystal 



field, the six spinful t2g orbitals break into a doublet 
|F+,±^) = -^{-iyz\Ut)^zx\i,t)Tixy\tA)), and 
a quartet |r+,±^) = {z^iyz - zx) \ i,^), |r+,±f) = 
^ {±iyz\ t,;) + zx\ t,;) + 2ixy\ ;,t))[32]. 

lifts their degeneracy at the F point and is diagonal 
in the F^ ,Fg" basis with eigenvalues of ^ and — ^/2, re- 
spectively. We set ^ = 2Ao/3 = 19.3 meV, which leads 
to the same orbital splitting as the DFT results. In the 
original t2g basis {yz\ t), yz\ i), zx\t), zx\ I), xy\t), 
xy\D), has off-diagonal terms and reads 
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+ is the TB Hamiltonian of bulk STO including 
SOC. Its bandstructure is shown as a solid line in Fig. 1(b) 
and agrees well with the DFT results. This model al- 
lows for a deeper understanding of the SOC effects: the 
SOC eigenstates are admixtures of the yz^ zx and xy 
orbitals, which explains the significant changes of the ef- 
fective masses. 

DFT results for interfaces: The band structure of 
LAO/STO (1.5/6.5) calculated by DFT is shown in 
Fig.2(a). Without SOC, similar to bulk STO, ah bands 
exhibit a parabolic-like behavior. In the x direction, the 
yz band is the fiattest (heaviest); at F, it is degener- 
ate with the zx. Due to the interface, the xy band is 
A/ = 250 meV lower in energy at F than the degenerate 
yz and zx bands. The splitting A/ is the most notable 
feature of the heterostructure[22, 23, 26]. It is not mainly 
a crystal field effect, but originates from the vanishing of 
the hopping from the interface Ti yz to LAO along the 
z direction [31]. Consistently, a similar behavior is ex- 
pected in a SrO terminated STO surface, and indeed we 
obtain it with A/ = 320 meV, see Fig. 2(b). Our cal- 
culated splittings are qualitatively consistent with the 
ARPES measurements of STO surfaces [33, 34]. 

Including SOC does not influence the xy band very 
much. It splits the degenerate yz and zx orbitals with 
Ao = 19 meV at F, see Fig. 2(a). The ab initio calculated 
Ao is qualitatively consistent with experiment, albeit 
smaller than its experimental value [33] Ao = 60 meV. 
Around F, the effective masses for xy and the resulting 
two states are 0.48, 1.14 and 0.72 mg, respectively. For 
the asymmetric case Fig. 2(c) the SOC also results in a 
spin splitting which is most noticeable at the xy-yz cross- 
ing region where it is up to 18 meV, see Fig.2(d,e). This 
spin splitting is a multi-orbital effect, very different from 
the standard Rashba spin splitting of single orbital. For 
a better understanding, we now construct a TB Hamil- 
tonian. 

Spin splitting at the interface layer: Without SOC, a 
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FIG. 2: DFT calculated band structure of (a) LAO/STO with 
two symmetric n-type interfaces, (b) SrO terminated STO 
surface, and (c) vacuum/LAO/STO with a single n-type in- 
terface. The region with biggest spin splitting around the 
xy-yz crossing point is magnified in (d) and for a LAO/STO 
(4/4) with asymmetric n-p interfaces in (e). A/ denotes the 
interface induced orbital splitting energy; the black dashed 
line is without SOC; the red and blue solid lines are the two 
opposite spin channels with SOC. 



model Hamiltonian H'q can describe the interface hopping 
and the induced splitting A/. In contrast to H^^ the hop- 
ping terms of in direction z essentially vanish. The 
diagonal term for xy is hence e^^ — 2ti cos kx — 2ti cos ky — 
^2—4^3 cos kx cos ky^ while that for the yz (and zx) orbital 
is e^^ — 2t2 cos kx — 2ti cos ky —ti — 2ts cos ky. The local 
energy terms e^y^y^ will be influenced by the interface 
crystal field, electron filling and confinement [23, 26, 31]. 
For simplicity, we approximate these by the bulk value 
Sq. Thus, A/ = ti — ^2 + 2^3 =0.4eV which is comparable 
to the DFT results. At the interface the Oh symmetry 
breaks down to C4^, and we can use the same atomic SOC 
matrix as before, since under the F^ doublet does 
not break whereas the Fg" quartet breaks into Fg" F^ , 
with the same set of basis functions as given previously. 
The Hq + Hamiltonian gives an atomic SOC induced 
orbital splitting about Aq = at F similar to Figs. 

2 (a-c). However Hq + does not contain any terms 
breaking the interface inversion symmetry, and hence it 
does not include the Rashba spin splitting. 

To this end, we introduce a term to describe the 
broken inversion symmetry at the interface, a key com- 
ponent for Rashba spin splitting. The essential physics 
of this term was analyzed by Lashell et al. [16] and then 
introduced by Petersen et al. [35] to construct a TB model 
for the Rashba effects of s-p orbit als in metal surfaces. 
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FIG. 3: (a) Schematic figure of the orbital deformation at 
the interface resulting in an anti- symmetric hopping term 7. 
Band structures calculated by (b) TB Hamiltonian Hq-\-H^-\- 
H^, and by (c) DFT for comparison; (d) the standard Rashba 
spin splitting Ar can be observed after zooming in the band 
structure of panel (b); red and blue indicate opposite spin, 
(e) Spin-splitting in the lowest band of the TB Hamiltonian 
for different A/. The inset shows a A;a;-linear (-cubic) spin 
splitting for positive (negative) A/, (f) Angular dependence 
of the spin splitting at the standard Rashba region (green) 
and the xy-yz /zx crossing region (brown) at the k points 
displayed in (b). 



To our knowledge, there and in other publications, 
was always treated as a parameter and hence its utility 
and importance are strongly Hmited. In this study, we 
project the DFT results above onto maximally localized 
Wannier orbitals [31] and then directly extract the spin 
independent hopping term 

2i sin kx 

^' ' ' 

iky 



— 2i sin kx —2i sin ky 



describing inter-orbital hopping terms due to the in- 
terface asymmetry. The key hopping term is 7 = 
{xy\H\yz{R)), where R is the nearest neighbor in x di- 
rection. As shown in the schematic Fig. 3 (a), 7 is an anti- 
symmetric hopping between xy and yz orbitals along the 
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X direction. Its origin is the interface asymmetry de- 
forming the orbital lobes of the interface layer. We find 
7 ~ 20meV at the n-type interface for all geometries, and 
hence take this value in the model. Let us note 7 drops 
quickly in the second and further layers towards its bulk 
value 7 = 0. 

The combined model Hamiltonian -\- -\- H^^ in- 
cluding Rashba effects, is expressed in the t2g basis by 
a 6 X 6 matrix, where Hq describes the interface hop- 
ping and splitting A/, includes the atomic SOC 
and accounts for the orbital splitting of Aq, and 
describes the interface asymmetry. The first effect is 
a standard Rashba-type of spin splitting in the single 
xy band. It splits a single parabola around the min- 
imum at r into two parabolas with opposite spin, see 
Fig. 3(d). By downfolding the matrix onto an effective 
Hamiltonian for the xy band, we obtain an analyti- 
cal expression for the spin splitting A^^ = 2aRkx with 
an = 2a^7/A/ = 0.76 x IQ-^eVA for A/ = 0.4eV, 
^ = 19.3meV, 7 = 20meV. Note that A/ depends 
strongly on the details of the interface and hence aR can 
be up to 8 times larger at A/ = where the formula 
above does not hold anymore, see Fig. 3(e). This well 
agrees with the experimental magnitude of the Rashba 
spin splitting = 1 — 5 x 10~^eVA[l]. 

If we turn A/ negative, which is possible by inter- 
face engineering [31], the lowest band is a mixture of 
yz and zx. In this situation there is no standard /c-hnear 
Rashba effect any more but we obtain a spin splitting 
2askl with 0^3 = 4eVA^, see Fig.3(e). Hence, the TB 
model also explains qualitatively and quantitatively the 
unusual /c-cubic spin splitting reported in [36] in a single 
framework, reconciling the puzzling discrepancy between 
experiments [36] and [1]. 

An even much bigger spin splitting 18meV occurs how- 
ever at the xy-yz crossing point, see Figs.3(b,c,e) and 
Figs.2(c-e). Unlike the isotropic splitting A^^ around 
F, this spin splitting is not only much larger but also 
anisotropic, see Fig. 3(f). This multiband effect is a par- 
ticularity of transition metal oxides and not occurring 
in semiconductors or metal surfaces. Experimentally, a 
similar anisotropy has been observed in LAO/STO het- 
erostructures with a particular strong SOC effect [37]. 

An important aspects of our study is also that the ex- 
ternal electric field [1, 2] does not significantly tune the 
SOC directly. As mentioned in the introduction, its di- 
rect contribution is too small. Even without it, we obtain 
the correct magnitude of the spin splitting. Nonethe- 
less, the spin splitting depends on the electric field [1, 2]. 
The explanation for this is an indirect effect: the elec- 
tric field tunes the carrier densities [38-42], band filling 
[23, 26, 43, 44], and the effective 7. The multi-orbital 
complexity might account for the discrepancy of the two 
reported spin energies tuned by gate voltages [1, 2]. 

Conclusion: We performed first principle calculations 
and developed a realistic three-band {xy^ yz and zx) 



model for SOC effects at LAO/STO interfaces and STO 
surfaces. The key ingredients to the spin splitting are the 
atomic SOC and the interface asymmetry, which enters 
via asymmetric t2g orbital lobes. The xy orbital around 
F exhibits the standard Rashba spin splitting 2aRk with 
aji = 2a^7/A/ ~ 10~^ eVA; in contrast, for negative 
A/ there is instead a /c-cubic dependence spin splitting 
in the lowest band around F. As A/ depends on the par- 
ticular surface or interface, this solves the experimental 
controversy regarding linear or cubic Rashba splitting. 
Even more importantly, we find an unusually large spin 
splitting 18meV at the crossing point of xy and yz/zx 
orbitals. Our results indicate that LAO/STO has pecu- 
liar SOC properties arising from the multi-orbit als char- 
acter which are absent in the standard single-band de- 
scription as for the nearly free 2DEG in semiconductor 
heterostructures. 
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